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ABSTRACT 

'66 ; 

■ 

■ We present the first results from tfie analysis of GIRAFFE spectra of more than 
^ . 1200 red giants stars in 19 Galactic Globular Clusters (GCs), to study the chemical 

^ i composition of second generation stars and their link with global cluster parameters. 

I We confirm that the extension of the Na-0 anticorrelation (the most striking signature 

I of polluted, second generation populations) is strictly related to the very blue (and hot) 

p\J ' extreme of the Horizontal Branch (HE). Long anticorrelations seem to require large mass 

■ and large-sized, eccentric orbits, taking the GCs far away from the central regions of the 

■ Galaxy. We can separate three populations in each cluster (primordial, intermediate 
I . and extreme) based on the chemical composition. In all GCs we observe a population of 

2 . primordial composition, similar to field stars of similar metallicity. We find that in all 

c/3 I GCs the bulk (from 50 to 70%) of stars belong to the intermediate component. Finally, 

1 the extreme, very oxygen-poor component is observed preferentially in massive clusters, 
but is not present in all massive GCs. 

> : 

■ Subject headings: stars: abundances — stars: evolution — globular clusters: general 

vn ■ 
rn ; 

. 1. Introduction 

oo ; 

■ Galactic Globular Clusters (GCs) show a series of star-to-star abundance variations for several 
^ ' light elements. Part of them (Li, C, N) show the same behavior in field stars, some others (O, 

^ . Na, Mg, Al) show a behavior peculiar to GCs. In particular, there are anticorrelations between 

O and Na, Al and Mg: they are present in all evolutionary phases sampled (both unevolved main 
sequence and evolved red giant branch stars) in all GCs studied so far (for a recent review see 
Gratton, Sneden & Carretta 2004). 

The fact that these anticorrelations are present also in unevolved stars (Gratton et al. 2001, 
Ramirez &: Cohen 2002, Carretta et al. 2004), unable to activate the required nucleosynthetic 
chains (p-capture in H burning at very high temperatures, e.g. Denisenkov & Denisenkova 1989) 
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implies that they have been imprinted in the gas by a previous generation and that we are seing 
a second generation of stars in these GCs. We do not yet know what kind of stars produced the 
pollution; the most popular candidates are intermediate-mass AGB stars (e.g. D'Antona & Ventura 
2007), or very massive, rotating stars (e.g.Decressin et al. 2007). 

The Na-0 anticorrelation is also stricly connected to the helium abundance (Na-rich, 0-poor 
stars should also be He-rich) and to the multiple sequences that have been recently found in some 
GCs (w Cen, NGC 2808, see Bedin et al. 2004 and Piotto et al. 2007, respectively) and attributed 
to populations with different Y. 

There are currently data on more than 1400 stars, in about 20% of known GCs, thanks to 
the efforts of many researchers, in primis the Lick- Texas group, lead by Chris Sneden and Robert 
Kraft; our group is now one of the main contributors, with more than 1200 stars analyzed. The 
bottom line is that we can stop talking about "abundance anomalies" and we must start speaking 
about the normal chemical composition of second generation stars in GCs. Times are ripe to start 
a more quantitative approach. 

Therefore, with the goal of understanding what is the relation between these abundance pat- 
terns and the horizontal branch (HE) morphology or structural parameters of the GCs, we have 
conducted an homogeneous survey of RGB stars in 150 Galactic GCs using the multiobject spec- 
trograph FLAMES@VLT. We obtained Giraffe spectra (at R~20000, comprising the Na I 568.2- 
568. 8nm, 615.4-6.0nm and [O I] 630nm lines) of about 100 stars per cluster, and UVES spectra (at 
R~40000, covering the 480-680nm region) of about 10 to 14 stars per cluster. 

We have targeted GCs with [Fe/H] ranging from —2.4 to —0.4 dex, with very different HB 
morphology (stubby red HB, blue HB, very extendend blue HB, bimodal HB), and very different 
global parameters (total mass, concentration, age, orbit, etc.). 

We have determined atmospheric parameters and measured Fe, O, and Na abundances for 
more than 2000 stars (of which more than 1200 cluster members and with both O and Na detected. 

In Fig. 1 we show a collage of the Na-0 anticorrelations observed in all 19 clusters in our sample. 
Solid lines separate the Primordial, Intermediate and Extreme populations, discussed below. 

The number of stars measured in each cluster depends on the richness of population, metallicity, 
S/N ratio and in some cases on field stars contamination, as in the case of the bulge clusters 
NGC 6388 and NGC 6441, or the disk cluster NGC 6171. 



1 NGC 104, NGC 288, NGC 1904, NGC 2808, NGC 3201, NGC 4590, NGC 5904, NGC 6121, NGC 6171, 
NGC 6218, NGC 6254, NGC 6388, NGC 6397, NGC 6441, NGC 6752, NGC 6809, NGC 6838, NGC 7078, and 
NGC 7099 
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2. Extension of the Na-O anticorrelation and (blue) HB 

As a way to quantify the extension of the anticorrelation, we adopted the interquartile range 
(the middle 50% of the data) of the [0/Na] distribution. This was suggested by Carretta (2006) as 
the best indicator, because it is not much sensitive to the ouliers in the distribution. 

The most tight correlation we found (shown in Fig. 2) is between the extension of the Na-0 
anticorrelation (well measured by the Inter-Quartile range IQR[0/Na]) and of the HB (measured 
by the maximum Teff reached on the HB). This has been first demonstrated by Carretta et al. 
(2007d) on a subsample of our 19 GCs and literature data, and amply confirmed on our whole 
sample. In Fig. 2 the red dots are for our data, the blue triangles for literature data: M 15 (Sneden 
et al. 1997), NGC 362 (Shetrone & Keane 2000), and M 13 (Sneden et al. 2004, Cohen & Melendez 
2005). 

M 15 is somewhat uncertain: the value from literature falls very well on the relation defined 
by all other clusters, while our value is more discrepant: maybe we missed some very 0-poor stars, 
difficult to measure at this very low metallicity. The Spearman rank coefficient tells us that the 
probability to get such a tigh relation only by chance is neghgible, ranging from 1 in 10~^ if we 
include our value for M 15, down to 1 in 10""^ or 4 in 10~® excluding M 15 or taking the Lick value. 

We conclude that this is a real, very strong relation: the first conclusion we can draw is that a 
same mechanism drives/affects the extent of the pollution on the RGB and the morphology of the 
bluest end of the HB. 

Notice, however, that this can be only considered a second parameter: the global distribution 
of stars on the HB (as indicated e.g., by the HB ratio (B+R)/(B+R+V) -see Harris 1996 and web 
updates) is not correlated with the IQR[0/Na], i.e. with the extension of the Na-0 anticorrelation. 

3. Relations with structural and orbital parameters 

We also found a good correlation between the IQR[0/Na] and the total mass (or absolute 
magnitude My) of the GCs (see Fig. 3): an high mass seems to be a requisite for an extended 
anticorrelation. 

This may be explained by a better capability of the cluster to retain the gas polluted by the 
first-generation stars to form the second-generation stars in which we measure the anticorrelated Na 
and O. However, high mass is not a sufficient condition: 47 Tuc is a notable counter-example. This 
cluster simply does not show very 0-poor stars and presents a short/normal Na-0 anticorrelation 
even if it is a very massive object. Some other factors must be involved. On the other hand, we 
note that the position of GCs standing leftward of the main relation (namely NGC 288, NGC 6218 
and M 71) is more easily explainable by the fact that they may have lost a (probably large) fraction 
of their mass in the past. 
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We explored the relations with orbital parameters, as first suggested by Carretta (2006), and 

performed simple bivariate correlations, adding to the cluster absolute magnitude a dependence on 
the cluster orbital parameters. In Fig. 4 we show the IQR[0/Na] as a function of a combination of 
My and the eccentricity of the orbits, whenever available (Dinescu et al. 1999). 

When we combine mass (magnitude) and orbital parameters the relations get tighter and 
statistically more significant. More importantly, 47 Tuc does not stand out at all in this plot, and 
is now virtually indiscernible from the other clusters. The same holds using other orbital parameters 
(as total energy, or period of the orbit, or maximum height over the plane). 

We consider these results as an indication that also the initial position and the time spent 
away from possible perturbations affect the possibility of a GC to build a gas reservoir from which 
the second-generation stars are born. 

4. The fraction of stars with Primordial, Intermediate and Extreme composition 

The ratio of first to second-generation stars is very useful to put constraints on the formation 

scenario. We assume as first-generation stars (or primordial, P) those with O and Na similar to 
field stars of same mctallicity. Wc divide the second-generation stars according to the pollution 
degree in Intermediate (I) and Extreme (E: those with [0/Na]< —0.9 dex, from a comparison with 
the distribution in NGC 2808, see Carretta et al. 2006) 

Wc have applied this operative division to all the 19 GCs in our sample. The lines separating 
the three components are shown in Fig. 1. From this figure and computing the fraction of stars in 
each component we found that: 

• the Primordial population is present in all GCs. 

• the Intermediate 2nd generation constitutes the bulk (50-70%) of stars. 

• the Extreme 2nd generation is not present in all GCs. 

The fraction of Intermediate and Extreme second generation stars are anticorrelated, and the 
Intermediate population is anticorrelated with the IQR([0/Na] (Fig. 5) since, by definition, the 
extension of the Na-0 anticorrelation is driven by the Extreme population (Fig. 6). 

When the fractions of the three components are plotted vs total mass (measured by total My, 
Fig. 7), we see that less Primordial stars are lost by the more massive clusters, as expected. As 
with the IQR[0/Na], the Extreme population is more pronounced in heavier GCs, but the relation 
is not one-to-one. 

The relations (also shown in Fig. 7) with half- mass radius (r^) are very tight, indicating that 
the shape of the anticorrelation was inprinted in the early phases of the cluster life, since is one 
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of the parameters less sensitive to changes due to the cluster evolution, both internal and due to 
external sources. More detailed studies, involving multivariate analysis, are under way. 

At the moment, our preliminary findings can be summarized as follows: 

• the extent of the Na-0 anticorrelation is strictly related only to the very blue and hot extreme 
of the HB 

• long anticorrelations are seen in clusters with both high mass and favourable orbital param- 
eters, taking them for long periods away from the central regions of the Galaxy 

• in all GCs there is a surviving population of stars of primordial composition, but 

• the bulk of the stellar population in GCs is composed of stars with intermediate degrees of 
pollution 

• finally, an extremely 0-poor component is observed only in some clusters and at the moment 
is not entirely clear why. 

In the near future we will present the full analysis of the clusters not already publishecl^, a 
systematic study of the properties of second generation stars and their relations with global cluster 
parameters and the relation between the Na-0 anticorrelation and HB morphology (Carretta et al. 
2009a, b and Gratton et al. 2009, in preparation, respectively). 

This work is based on the data collected at ESO telescopes under the programmes 072.D-0507 
and 073.D-0211. We acknowledge partial financial support from PRIN-INAF 2005 "Experimenting 
stellar nucleosynthesis in clean environments" 
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Fig. 1. — Summary of the Na-0 anticorrelation observed in the 19 globular clusters of our sample. 
Arrows indicate upper limits in oxygen abundances. The 3 solid lines in each panel separate 
the Primordial component (located in the Na-poor/O-rich region), the Na-rich/O-poor Extreme 
component, and the Intermediate component in between (see text for more details). 
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Fig. 2. — Interquartile range of the [0/Na] ratio as a function of the maximum temperature reached 
on the HB (taken from Recio-Blanco et al. 2006). Fihed red circles are the GCs of our sample, 
filled blue triangles are estimates for 3 clusters from the literature. The solid line connects the 
values for M 15 (see text). 
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Fig. 3. — The extension of the Na-0 anticorrelation (measured using the IQR[0/Na]) as a function 
of the clusters' absolute magnitude My from Harris (1996 -as updated in the web page). The 
Pearson and Spearman correlation coefficients are listed. 
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Fig. 4. — The extension of the Na-0 anticorrelation (measured using the IQR[0/Na]) as a function 
of a combination of the clusters' absolute magnitude My and of the eccentricity of the orbits in the 
Galaxy. 
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Fig. 5. — The extension of the Na-0 anticorrelation as a function of the fraction of stars in the 
intermediate I component of second-generation stars. 
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Fig. 6. — The extension of the Na-0 anticorrelation as a function of the fraction of stars in the 
extreme E component of second-generation stars. 
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Fig. 7. — Left: the fraction of stars in the P, I and E components of second- generation stars as a 
function of the total absolute magnitude (hence, mass) of clusters. Right: the same, as a function 
of the half-mass radius of the clusters. 



